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THE  PROBLEM 


Provide  in  u  single  reference  the  analytical  expressions  for  the  response 
of  transducer  arrays  of  standard  form.  Establish  the  most  imiiortaut  factors  in  se¬ 
lecting  a  particular  configuration  for  target  detection  and  or  bearing  estimation. 


RESULTS 

1.  Analytical  expressions  of  responses  foi  linear,  planar,  cylindrical, 
and  spherical  arrays  are  given. 

2.  Additive  arrays  are  “best.”  for  detection  of  a  point  target  in  a  uniform 
noise  field. 

3.  Split-beam  multiplicative  arrays  are  “best”  for  bearing  estimates  of 
a  single  point,  target. 


RECOMMENDATION 

Use  the  analytical  expressions  of  array  response  as  a  basis  for  computer 
programs  used  in  predicting  sonar  performance. 
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INTRODUCTION 


In  the  continuing  effort  to  improve  the  capabilities  of  sonar  systems,  the 
most  common  form  of  signal  processing,  and  one  of  the  most  useful,  has  lx*on  by 
the  use1  of  directional  transmitting  and  receiving  transducer  arrays.  The  array 
arrangement  makes  it  possible  to  place  greater  [lower  in  the  water  and  direct  that 
[lower  as  desired.  Directivity  offers  important  advantages  over  omnidirectional 
systems.  In  the  transmitting  transducer,  directivity  provides  a  greater  percent  of 
the  output  energy  in  a  desired  direction  (in  the  main  beam)  and  thus  provides  more 
energy  on  target  (if  it  is  in  the  main  beam)  than  would  lie  provided  by  an  omnidi¬ 
rectional  transducer.  In  the  receiving  array,  the  directivity  discriminates  against 
any  noise  sources  which  are  not  in  the  main  beam,  and  thus  improves  the  signal- 
to-noise  ratio.  Considering  the  reciprocity  of  antenna  systems,  all  transducer 
systems  may  he  considered  to  operate  as  receivers. 

In  the  study  to  be  reported  here,  arrays  of  isotropic  elements  were  used 
to  provide  directionality  and  so  improve  detection  probability  and  parameter 
estimation  -  the  two  major  aspects  of  signal  processing.  For  detection,  the  pri¬ 
mary  problem  is  to  maximize  the  signal-to-noise  ratio  and  at  the  same  time  meet 
other  system  requirements  of  size,  frequency,  bandwidth,  scan  time,  time  on  target, 
etc.  For  parameter  estimation  (in  this  case  direction  of  arrival  of  the  wave  front), 
improved  sensitivity  of  the  output  to  changes  in  the  parameter  is  of  importance, 
as  well  as  signa’-to-noise  ratio. 

The  following  section  descrilies  various  arrays  of  isotropic  elements  to 
provide  directionality.  Subsequent  sections  are  concerned  with  shading  or  weight¬ 
ing  of  the  elements  to  reduce  side-Iolic  level,  and  with  some  of  the  common  forms 
of  tracking  systems  used  for  hearing  estimation. 


ARRAYS  OF  ISOTROPIC  ELEMENTS 


Linear  Arrays 

The  first  array  considered  is  a  line  of  equally  spaced  isotropic  elements 
(the  response  of  an  isotropic  element  is  constant  in  all  directions)  as  shown  in 
figure  1.  To  obtain  the  fur-field  response  due  to  a  plane  wave  arriving  at  the 
array  from  an  angle  0  relative  to  the  normal  to  the  array,  the  output  of  each  ele¬ 
ment  is  determined  and  then  summed: 

R(0)  -  y  exp  (i  t!/k) 
k 

where 

R(0)  is  the  array  response 

R k  is  the  amplitude  response  of  the  ft11'  element 
\bk  is  the  phase  at  the  element 

Since  only  the  far-field  response  is  being  considered,  the  amplitude  at  each  cle¬ 
ment  is  assumed  equal  and  the  response  becomes 


Referring  again  to  figure  1,  the  phase  at  the  kth  element  relative  to  the  center 
element  is  given  by 


=  d.  cos  6 = kdr  cos  6 
A  A 


where  d  is  the  element  spacing,  dr  =  - , 

A 


and  A  is  the  wavelength  of  the  source.  Thus 


R(d)  =  A  exp  {ikd rcos  6) 


The  sum  from  -k0  to  k0  is  of  the  form 


^  ^  exp  (/M)  =  exp  (-/&<£) 


+  exp  t- 1  (ft<>  - 1)  6  +  •  •  • +  1 +  •  •  •  +  exp  [/  (k0  - 1)6]  +  exp  Uk06) 


.exp(ikod>)  •  —  =  l*exp(i<6)  +  exp(/26)  +  -  •  ■  +  exp(/2/f0<£) 
A 


Therefore 


l-exp[/'(2/o0*  1)6] 

exp{ik06)  -  =  — ; - - - 

A  l-exp</6) 


/2k0+  1\  .  1/2ft0-l\ 

R  cx|)(-ife06)oxp[i(fc. ^  l)d]  LJ~ 

A"  1-«PW  '  oxp(-.|)-exp  (i|) 


Since  2/j0  +  1  =  n.  the  total  number  of  elements. 


R  exp(-d|)-cxp(.d|)  sin^ 
A  /  •  <5\  /-  <£  \  •  <5 

cxp\  T/  expvT/  sinT 
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and 


RUH 


.  \  s  >  n 


■) 


cos  0 


COS  t) 


if  .A  is  chosen  to  normalize  the  output  to  1  whoa  0  =  0,  then  A  -  I  n  and 


rid 

sin  — -  cos 0 

RW)  , - = - 

d  r 

1 1  sin  —  cos n 

If  a  phase  delay  on  is  added  to  each  element  output  such  that  the  delay  for  the 
/i1h  element  is  ho„  the  response  for  the  array  will  be  given  by 


RW) 


ln(lr  rcd>- 

sinl— —  cos  d  -  —j— 


■  (dr  a  \ 
n  sin^—  cos  0  -  — ) 


The  proper  phase  delay  between  adjacent  elements  *o  steer  the  center  of  the  main 
beam  to  d0  is  ci0  d,  cos  0  and  the  response  is  given  by 


.  ndr 

sin  — — (cos  0- cos  d0) 

R’O,  0n)  - = - 

dr 

n  sin —  (cos  li  -  cos  (L) 
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I  he  response  RUi,i)„)  is  an  amplitude  pattern  or  directivity  function  or  array  factor, 
in  determining  directivity  index  and  signal -to-uoise  power  at  the  receivei,  the 
power  pattern  :s  often  ot  interest.  The  power  pattern  is  proportional  to  the  square 
of  the  magnitude  of  the  directivity  function: 


PUD  h'  IrudU 


Planar  Arrays 


The  directivity  function  of  an  array  of  isotropic  elements  is  given  by: 


R  =  Rk  exp (6k) 


where  R ^  is  the  amplitude  of  the  fcth  element 

<t>k  is  the  phase  of  the  fcth  element  referenced  to  the  origin. 

The  phase  term  is  given  by 


,  2  n  , 

where  A  is  the  difference  in  the  distance  from  the  source  to  the  origin,  D0,  and 
the  distance  from  the  source  to  the  kth  element,  Dk: 

Da1  -  Dfe’  Do*  -  Dfc> 

A'D«-Dfc  =  0o +  =  2Dd  -  A 

If  the  source  has  coordinates  x,  y,  z  and  the  kth  element  x*,  yk,  zk  then 
x  =  D0  cos  ox  ,  xk  ~  Pk  cos  axk 

y  =  Do  cos  ay  yk  =  pk  cos  ayk 

z  =  D0  cos  a  z  zk  =  PkCOSazk 

with  pk  the  distance  from  the  origin  to  the  ftth  element  and  ox,  av,  a z  the  angles 
shown  in  figure  2.  Using  these  equations 


A  ^  2073a  ID°a '  <x ' Xfe)>  ■ (y ' y*)l  ■ (* "  ^)31 

=  2o^a  \2xxk  +  2yyk  +  2z2k- {v +  -V +  v>l 

2D0  pk  cos  ax  cos  axk  +  2 D0  ;>k  cos  ay  cos  ayk  +  2 D0  pk  cos  az  cos  azk  -  pk 1 
=  2D0- A 

Pk  \  Pk  1 

=  lcos  a* cos  a*k  +  cos  ay cos  ayk  +  cos  cos  a*k '  2Do\ 

2 Do 
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and  the  response  is  specified  by  the  product  of  the  response  in  the  Y-Z  plane 
and  the  response  in  the  Y-Z  plane.  To  use  this  separation  of  variables  it  is 
important  to  use  the  proper  coordinate  system.  As  an  example,  choose  the  0 ,  6 
system  shown  in  figure  3  and  transform  the  response: 


cos  ax  =  sin  <f> 
cos  av  =  cos  (f)  cos  d 


which 

single 


R(6 ,  <£)  = 


(Mnd  \  find  \ 

I— ^  sincS )  sin  I— —  cos  0  cos 

Hx  .  \ 

iin  1——  sin  cf>l 


(nd 

N  sin(-^cos<£cos0J 


is  no  longer  a  product  of  two  functions  each  of  which  is  a  function  of  a 
variable. 


x 


Figure  3.  Geometry  for  azimuth-elevation  coordinate  system. 


Introducing  a  phase  shift  between  elements  to  steer  the  main  beam  to  a 
direction  cosaJCS,  cos  ays  and  writing  the  response  in  terms  of  the  differences 
of  the  direction  cosines  results  in  the  following: 


12 


Figure  4  is  a  plot  of  these  equations  near  the  main  beam.  Note  that  for 
a  square  array  ( M  =  N )  the  main  beam  is  circular  in  cross  section  at  the  half-power 
points,  and  it  is  found  to  lx?  elliptical  for  ,\l  £  N.  As  R  approaches  the  first  null, 


a  broadening  of  ihe  Ix'nm  along  the  diagonal  is  noted.  The  half-power  beam  width 
of  a  uniformly  i  1  i uni inat 1  square  array  is  given  in  r-spaee  as 


C„  0.88H 

.1  yd  (length  of  the  aperture) 

Other  cixdTicienls  are  applicable  for  other  ty|X>s  of  ajxirture  illumination. 
In  general,  the  main  lx.'am  broadens  and  the  side-lobe  level  is  reduced  as  the 
illumination  function  is  tapered.  Further  discussion  of  tapering  or  shading  the 
illumination  will  appear  later  in  the  report. 

The  main  advantage  in  using  the  r -space  is  that  the  lx?um  shape  is  not  a 
function  of  the  scan  position  of  the  Ix-arn  when  [dotted  in  r -space.  The  lx;am 
maximum  is  scanned  by  introducing  a  phase  difference  i b  and  uv  IxHween  ele¬ 
ments  in  the  array.  The  direction  cosines  corresponding  to  the  phase  shifts  are 


If  a  complex  T -plane  is  defined  by 


T  COS  a  I  i  COS  rt , 


s  COSaxx  ’  '  C<>8  “vs 


1.7  O 

r  ~ ,  i  ii. 

•'  dr  dr  dr 


1.7  U  •  /  (.  / 


1  lie  complex  .--plane  now  become 


i  r_  7  T 


The  array  response,  which  is  a  function  of  the  complex  variable  r,  is  in¬ 
variant  in  the  T-plane.  Scanning  the  beam  by  introducing  the  complex  phase  delay 
•Ir  =  ifr x  ■*-  i  ilr  simply  translates  the  pattern  in  space  so  that  the  pattern  center 

moves  to  T  (See  fig.  5.)  If  the  phase  delay  is  chosen  such  that  ]T|>1  the 

s  a 

half-power  contour  moves  outside  the  unit  circle  and  becomes  imaginary  (i.e.,  un¬ 
observable). 


COS  a  v 


Figure  5.  Scanning  in  the  complex  T  plune  with  one 
wavelength  element  spacing. 


The  effects  of  scanning  the  main  Ijeain  may  now  lx?  considered  in  the  T  ■ 
plane  with  only  a  translation  of  the  pattern  being  considered.  Figure  5  shows  the 
main  !x?am  and  the  secondary  maxima  of  a  square  array  of  isotropic  elements.  The 

secondary  maxima  are  spaced  rT0  -  ,-\0  ~r  apart  and  there  are  an  infinite  numljor 

of  them.  It  must  lx?  remembered  that  only  those  contained  within  the  unit  circle 
are  real.  The  larger  the  element  space  the  closer  the  spacing  of  the  secondary 
maxima  and  hence  the  greater  numlx?r  contained  in  the  unit  circle.  If  the  elements 
are  spaced  one  wavelength  apart,  the  main  lx?am  and  four  end-fired  secondary 
beams  are  visible  and  any  scanning  from  broadside  will  move  the  end-fired  beams 

within  the  unit  circle.  Closer  spacing  of  the  elements  can  suppress  the  secondary 

2A 

beams.  For  example,  with  d  ^  no  secondary  maxima  occur  if  T s  <  1  2  (fig.  {->). 

Further  suppression  of  secondary  maxima  is  jrossiblc  when  using  nonisotropic 
sources  with  low  gain  in  the  direction  of  the  expected  secondary  maxima. 


COS  a  v 


Figure  6.  Secondary  maxima  of  a  scanned  array  with  two-thirds 
wavelength  element  spacing. 


The  value  of  using  the  T-plane  or  direction  cosine  space  when  examining 
the  effects  of  scanning  the  main  beam  of  the  array  has  been  shown.  However,  to 
permit  a  physical  interpretation,  a  transformation  to  a  spherical  coordinate  system 
is  necessary.  The  proper  transformation  is  found  by  expressing  the  complex  di¬ 
rection  cosine,  T,  in  polar  coordinates: 

T  =  cosax  +  »  cosay  =  sin  0  cos  6  +  i  sin  6  sin 

T  =  sin  9  (cos  <f>  +  i  sin  <£)  =  sin  6  exp (i'«5) 

Hence,  the  desired  transformation  is  the  projection  of  the  unit  circle  in  the  T- 
plane  onto  a  unit  sphere  (fig.  7).  By  limiting  the  values  of  0  to  0  <  6  <  £  only 

the  positive  hemisphere  is  of  interest.  In  the  case  of  a  rectangular  array  (M  i  N), 
the  transformation  is  the  same  but  the  main  and  secondary  maxima  will  be  approx¬ 
imately  elliptical  in  cross  section  rather  than  circular  as  discussed  earlier. 

Referring  to  figure  7,  as  the  main  beam  is  scanned  away  from  the  broad¬ 
side  position,  the  beam  width  changes.  The  reference  beam  width  at  broadside  is 

Bo  =  2  sin*1  (Ar) 


where 


Cl  Ao 


Thus,  as  the  beam  is  scanned  away  from  the  normal,  it  broadens  in  the 
0  direction  and  the  half-power  contours  for  a  square  array  change  from  circular  to 
elliptical. 

Because  of  the  change  in  the  beam  shape,  its  maximum  point  is  not  cen¬ 
tered  between  the  half-power  points.  A  beam  eccentricity  can  be  defined  by 

(ft -*.)-(*. -ft) 

"(ft -».)  +  (*. -ft) 

Using  the  above  equations  as  they  apply  to  the  half  beams  03  -  0S  and 
0S  -  0i  yields 

2e  =  tan  1/4  (0,  +  0,  +  20s)  tan  1/4  (0,  -  0t) 

Be 

Approximating  0!  +  0:=;2  0C  and  tan  1/4  (03  -  0,)  =  -— 

s  4 

yields 

Be 

e  -  —  tan  0„ 

8  s 

The  difference  between  the  beam  maximum  and  the  center  of  the  half-power  contour, 

ft..+  ft  ,  is  then 
2 

Be1 

A0  =  —  tan  0S 

Thus,  for  narrow  beams  (i.e.,  Bg  small)  the  eccentricity  may  be  neglected.  Figure 
8  is  a  plot  of  the  pointing  error  vs  scanning  angle  as  a  function  of  the  broadside 
beam  width. 


Figure  8.  Pointing  error  vs  senn  angle  ns  a  function 
of  the  broadside  t*ciim  width. 


Nonuniform  Illumination 

The  principles  which  have  been  developed  will  Ik?  applied  to  tapered  il¬ 
lumination.  The  side-lobe  level  of  a  uniformly  illuminated  array  is  13.2  dB  below 
the  main-beam  intensity.  To  achieve  lower  side-lobe  levels.  tapering  of  the  illum¬ 
ination  across  the  array  or  weighting  the  output  of  each  clement  is  used.  Many 
types  of  illumination  ta|x?rs  have  been  proposed  for  giving  various  relationships 
between  gain,  beam  width,  and  side-lobe  level.  The  choice  of  illumination  taper 
will  depend  upon  the  relative  emphasis  placed  on  these  three  parameters. 

The  uniform  distribution  gives  the  13.2-dB  side-lobe  level,  theoretically 
the  highest  array  gain  for  a  uniform  phase  distribution  and  moderate  lx?am  width. 

A  second  important  illumination  taper  is  shown  as  the  Dolph-Chebyshev 
distribution.  Using  Chebyshev  polynomials.  Dolph  has  synthesized  an  illumination 
distribution  giving  uniform  side  lobes  of  any  desired  amplitude.1  It  was  shown  by 
Dolph  that  this  distribution  gives  the  narrowest  beam  width  for  any  specified  side- 


'Superscript  numbers  identify  references  listed  at  end  of  report. 


I ( ) ' x '  level  and  conversely  gives  the  lowest  side  lobe  level  for  a  given  Iteaui  width. 
In  this  respect,  it  is  opt  iiiniin.  'The  distribution,  Ixdng  specified  by  polynomials, 
is  given  at  discrete  intervals  only  and  is,  therefore,  ideally  applicable  only  to 
arrays.  The  distrilmt  ion  is  characterized  by  current  peaks  in  the  end  elements 
which  Ix'come  inconveniently  high  for  very  long  arrays  with  very  low  side  lobes. 
For  this  reason,  the  |)olph-(  'hebysliev  distribution  is  usually  used  only  with  arrays 
of  moderate  length. 

Van  d('r  Maas'  adapted  the  Doiph-t 'hebysliev  distribution  to  a  continuous 
aportur*'  bv  allowing  the  number  of  elements  to  approach  infinity.  In  the  limit  the 
far-field  distribution  approaches  the  function 

F  (U,A)  -  cos  n 

where 


V 


Sd 

Ao 


cos  a  -  cos  a 


and 


17  =  cosh  n  A  is  the  side-lobe  ratio 


This  radiation  pattern  is  physically  unrealizable  because  of  the  requirement  for 
infinite  current  peaks  at  the  edges  of  the  aperature. 

Taylor1  has  modified  Van  der  Maas’s  distribution,  giving  slightly  greater 
beam  width  and  higher  gain  for  the  same  lobe  level.  The  far-field  pattern  for  a 
line  array  is  given  by4 


F(U,A,rj) 


Vr*  , _ u1 

r  (7;  •  I  )  1  ’  <77  —  (  ■’)  I  o'l  A3  +  (r/~  ‘/j)J] 


where 


7 _ 

[A! 


The  constant  7/  is  arbitrary  but  if  it  is  set  equal  to  cosh  n  A  then  F  <(./,  A,  17)  may 

Ixi  quite  closely  approximated  in  the  region  j— |  •  A  by  the  function 

a 


F  ( U,A )  =  cos  77  .4 ' - 1— j  1 

Using  this  approximation,  the  lieani  width  of  the  modified  Taylor  distribution  is 
greater  than  that  of  the  Dolph-Chebyshev  limit  by  the  factor  a.  Figure  V)  is  a  plot 
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A  number  of  other  illumination  tapers  are  frequently  used  when  moderate 
side  lobes  are  required  and  the  primary  requirement  is  convenience  and  ease  of 
calculation.  A  number  of  these,  such  as  “truncated  gaussian”  and  “cosine  squared 
on  a  pedestal,”  give  excellent  results  with  side  lobes  approximately  25  dB  down. 

A  comparison  of  the  change  in  beam  width,  pointing  error,  and  gain  will  now  be 
made  for  a  34  A  linear  array  at  scan  angles  of  30  to  55  degrees,  with  uniform  illu¬ 
mination  versus  a  modified  Taylor  distribution.  A  34  A  linear  array  of  isotropic 
elements  has  a  broadside  half-power  beam  width  of5 


e 


0 


0.888 

L/Ao 


=  0,0261  radian  =  1.496  degrees 


The  gain  at  broadside  is  * 


C  =  —  =  68  =  18.32  dB 
•  A® 

and  the  first  side  lobes  are  13.2  dB  below  the  main  beam.  When  scanned  to  30 
degrees  from  the  normal,  the  half-power  beam  width  is  increased  to 5 


cos  Vi  <0,  +  0,) 


where 


sin  0,  =  sin  30  -  sin 


•> 


sin  0,  =  sin  30  *  sin 
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•» 


Therefore 

0,  -  0.49372  radian 

0,  •-*  0.55398  radian 

0,o  -  0.0301427  radian  =  1.727  degrees 

The  point  of  maximum  gain  is  displaced  from  halfway  between  the  half-power 
points  by 


A  0,o 


0,oa 
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tan  30°=  6.55769  x  10"5  radians 


=  13.526  seconds 


and  the  gain  is  reduced  by  the  “projected  area” 

G,0  =  G  cos  0  =  58.89  =  17.70  dB 

'I  S 

or  a  loss  of  0.62  dB  compared  with  the  gain  at  broadside.  All  calculations  are 
based  on  isotropic  elements.  The  directivity  of  the  elements  may  be  taken  into 
account  as  outlined  by  Von  Aulock.5 

With  a  scan  angle  of  55  degrees  from  the  normal,  the  half-power  beam 

width  is 


cos  Vi  <#j  +  0,) 

6i  =  0.93764  radian 
=  0.98311  radian 

0, ,  =  0.0455  radian  =  2.607  degrees 
The  maximum  gain  is  displaced  by 

A0,,  -  -~-tan  55 
o 

A  0, ,  =  3.6959  x  10**  radians  =  0.021176  degree 

*  76.73  seconds 


The  gain  at  the  can  angle  of  55  degrees  is 

G,,  =  Gy  cos  ds  =  39.00  or  15.91  dB 

for  a  2.41-dB  reduction  compared  to  broadside. 

The  modified  Taylor  distribution  chosen  for  comparison  has  a  side-lobe 
level  of  30-dB  (power)  below  the  main  beam  and  rj  =  6.  The  ratio  of  the  beam 
width  of  the  main  lobe  compared  to  that  of  the  Dolph-Chebyshev  is  given  by 


[AJ  +  (^-Vi)J]  V, 


where 


rj  =  cosh  nA  is  the  voltage  side-lobe  ratio 


Thus. 


a  =  - 2 -  =  =  1 .0608 

[1.74229  -  (6-  lt)J]  ’j  5.656 

and  the  main  lobe  is  approximately  6  percent  wider  than  a  Dolph-Chcbvshev  (DO 
with  the  30-dB  side-lobe  level  which  is  given  by 


<)  =  60.6  —  (degrees) 

DC  L 


compared  to  that  of  a  uniform  illumination  given  by 

0  50.9  y  (degrees) 

Lj 

The  ratio  of  the  main  l>eam  of  a  modified  Taylor  distribution  with  side  lobes  30  dB 
down  to  a  uniform  illumination  with  side  lobes  13.2  dB  down  is 

P  =  •  SI  "  «  1.19056  =  1.262946 

50.9 

given  a  beam  width  at  broad  side  for  a  34  A  linear  array  of 
<)T  =  a  60.6  y  ( -30  dB  side  lobes) 

Lj 

=  1.890  degrees  =  0.03299  radian 

The  gain  of  the  modified  Taylor  distribution  relative  to  the  uniform  is 
given  by 7 

AC  2 

2  L  ~  1 .93  A  r f  *  2  Oj-1) 

for 

7)2  =  103 

A  =  1.32 
1  = 6 


yf  -  0.782  or -1.06  dB 

m  Lj 


or 


G  =  53.18  or  17.26 

At  a  scan  angle  of  30  degrees,  the  half-power  beam  width  has  increased  to 


Oio  =  2.18  degrees  =  0.038048  radian 


and  the  gain  decreased  to 


C30  -  46.05  or  16.63  dB 


The  displacement  of  the  maximum  gain  point  from  midway  between  the  half-power 
points  is 

0.  o2 

A  0ia  =  — —  tan  30  (radians) 


0.38038) 3 


(0.5774) 


=  1.0448  x  10'*  radians  =  21.55  seconds 

With  a  phase  delay  between  elements  to  give  a  55-degree  scan  angle,  the  half- 
power  beam  width  is  increased  to 

0t  =  0.9318 

6,  =  0.9893 

0, ,  -=  0.57519  radian  =  3.295  degrees 
The  gain  at  55-dcgree  scan  angle  is  reduced  to 

G, ,  =  53.18  cos  55  =  30.5  or  14.84  dB 
The  displacement  of  the  maximum  gain  point  is 


A0,,  =  x  1.4282 

8 


5.906  x  10'*  radians  -  121.83  seconds 


26 


Table  1  summarizes  the  above  results  for  a  linear  array  of  isotropic  elements  34 
wavelengths  long. 


TABLE  1.  SIDE-LOBE  LEVEL  FOR  UNIFORM  AND 
MODIFIED  TAYLOR  ILLUMINATION 


Uniform  Illumination: 

Side-lobe  level  (dB  down  from  main  beam) 

Half-power  beam  width  (degrees) 

(radians) 

Gain  (ratio) 

(dB) 

Displacement  of  beam  center  (radians) 

(arc  sec.) 

Modified  Taylor  Illumination  ( rj  -  6): 

Side-lobe  level  (dB  down  from  main  beam) 

Half-power  beam  width  (degrees) 

(radians) 

Gain  (ratio) 

(dB) 

Displacement  of  beam  center  (radians) 

(arc  sec.) 


Scan  Angle  from  Broadside  (c) 

0 

30 

55 

13.2 

13.2 

13.2 

1 .490 

1.727 

2.007 

0.0201 

0.03014 

0.0455 

08.0 

58.9 

39.0 

18.32 

17.70 

15.91 

0.558  y  H)  5 

30.959  >.  10“ 

13.53 

70.23 

30.0 

30.0 

30.0 

1.890 

2.18 

3.290 

0,0.3299 

0.03805 

0.057519 

53.18 

40.05 

•30.5 

17.20 

10.0.3 

14.84 

10.448  x  10  5 

59.00  ■<  10“ 

— 

21 .55 

121.83 

Cylindrical  Arrays 

Referring  to  figure  11.  the  response  or  directivity  function  for  a  cylindri¬ 
cal  array  is  given  by 


/?=  Y  iRJexp  2id>s[ 


where  l  is  an  index  for  the  stave  of  the  element  and  s  is  an  index  for  the  element 
location  on  the  stave  with  the  (0.0)  element  being  located  on  the  Z  axis. 

The  phase,  <£s{,  referenced  to  the  origin  is  given  by 

2 

*st  =-f  (xst  C0Sax+  >'st  C0Say  +  Zsl  COSaz) 


Figure  11.  Geometry  for  cylindrical  array.' 


The  coordinates  for  the  s,  t  element  are  given  by 

x  =  sds 

■  ltd‘) 

y  =  r  sin  yJT"  J 
z  -  r  cos 

dg  is  the  element  separation  in  the  x  direction 
dt  is  the  arc  distance  between  staves 


r  is  the  radius  of  the  cylinder 


If  phase  compensation  is  added  to  give  a  maximum  response  in  the  ( axs ,  a ys.  azs) 
direction  the  phase  of  the  (s,t)  element  is  given  by 


=  IT  [~T  <COSa*~  COSaxS)  +  sin  (~r)(coSay- 


COSavs\  +  cos 


(?) 


(COS  a  -  cos«?  \ 


For  the  special  case  of  compensation  for  a  plane  wave  along  the  Z-axis  the  phase  is 


The  first  term  in  the  equation  for  phase  with  general  compensation  is  the  only  term 
which  depends  on  the  index  s  and  it  is  in  the  same  form  as  the  response  for  a  lin¬ 
ear  array.  Thus  if  the  amplitude  term  is  not  a  function  of  s  the  response  may  be 
written  in  the  form 


R=Sexp("4*)'  2  lR(i  exp(/d>( 


■Rt  =  n  Rst' 


2n 


6s=~rSds  (COSa.v-cos«.rs) 


‘?„r  [  .  ldt  ,(^t  1 

6,  =  — —  Js  in  —  (cosav.  -  cos  ays)  +  cos  —  (C0Sa,-C0Sa.sJ 


and 


•  fnnd«  1 

D  V  S,nl  A  (C0S«X-C0S«.V,)J 

Rs-  )  exp(i6s)= - — - — J 

*  n  sin  —^-(cos^-cosU 


where  n  is  the  number  of  elements  in  a  stave.  The  response  of  the  cylindrical 
array  with  no  shading  along  the  axis  of  the  cylinder  may  be  viewed  as  the  product 
of  the  response  due  to  a  line  array  and  an  arc  array.  The  response  for  the  arc  of 
a  circle  is 


H 


^  nt  ,'XI,(,'S) 


;  ,7/ 


sin —  ms«vs)  cos —  (cos  «  ,  -  (  os  s\ 


For  the  special  case  where  the  array  is  compensated  for  a  plane  wave 
arriving  in  the  .v-g  plane  hut  with  a  tilt  angle  of  6r  the  response  is  given  by 

K  K,  •  K\ 

. ,  ,  1 


sin  — 

L  A 

•—  (cos,/v  sin6() 

n  sin  I  -• 

ul  . 

— i(co sin  6  . 

L 

A 

<‘M>(  <<!>,) 

Since  compensation  in  azimuth  is  normally  done  by  selecting  the  stave  in  the  de¬ 
sired  direction  as  the  center  of  the  are  the  above  ease  is  a  very  general  one. 
Figure  [2  is  an  example  of  a  pattern  calculated  for  the  ease  of 

<> 

» > 

“v  "  “ 
it  _  in  a  steps 


A 


lti  staves  separated  hv  7.7  around  tho  cylinder: 

S  omnidireetional  elements  per  stave  (s„  1; 

V’l  ■  1)  d, 

H,  eos  2 - - - 

1  2r 


The  plot  is  of  10  log  H  (,/.())  2  -  10  log  l\  (0.0)  Note  the  relatively 
large  response  to  the  rear  of  the  arrav.  This  mas  Ia‘  reduced  by  using  elements 
which  are  themselves  directional,  has  iug  sees  low  response  in  the  car  hemisphere. 


Figure  12.  Response  in  Y-Z  plane  for  cylindrical  array.. 


Spherical  Arrays 

The  equations  for  the  response  of  an  array  of  point  sources  is  easily 
extended  to  the  case  where  the  sources  are  arranged  on  the  surface  of  a  trans¬ 
parent  sphere  with  diameter  p.  The  general  equation  for  the  phase  of  the  frth 
element  with  direction  cosines  axk .  axk.  azk  and  distance  from  the  origin  pk  is 
given  by 


2  *  Pk 

6k  =  — —  (COSar/,COSaxJ-COSavfcCOSav  J-C0Sa;A  COSa;) 

where  the  reference  point  is  the  origin  and  the  incoming  plane  wave  has  direction 
cosines  ax.  av.  nz( see  fig.  13).  If  the  array  is  compensated  to  give  a  maximum 
response  in  the  (a  .  avs.  azs)  direction  the  resulting  phase  for  the  /elh  element 
is  given  by: 
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■Vs  I*'08".,/,  (COSf'.v  c’os«vS)  -  eo8«a(cos«v-cos  avs) 


►  COS  «  ,  A,  (  CDS  ,  -  COS  ,  s  )  J 


where  the  subscript  was  dropped  from  p  since  the  distance  from  the  origin  is  the 
same  for  each  element  of  the  spherical  array.  Because  of  the  large  number  of  pos¬ 
sible  arrangements  of  the  elements  on  the  sphere  and  the  fact  that  most  of  those 
arrangements  used  in  practice  do  not  result  in  any  great  simplification  of  the  equa¬ 
tions.  they  shall  be  left  in  the  general  form  as  given.  The  resulting  response  is 
given  by 


R  ^  Rf,[  exp (i6f.) 
k 

The  amplitude  term  is  used  to  account  for  the  sphere  not  being  transparent, 
the  directional  response  of  the  individual  elements,  and  any  shading  function  which 
is  incorporated. 
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ANGULAR  TRACKING  SYSTEMS 


Tlii*  characteristics  of  several  sonar  tracking  techniques  wort*  studied  lo 
determine  (ho  most  suitable  for  an  ASW  system  employing  a  planar  transdueer 
array.  The  choice  of  the  "Ix'st"  technique  from  those  available  was  approac  hed 
by  comparing  their  performance  trade-offs  in  applications  to  an  ASW  planar  array 
with  their  characteristics  in  applications  to  an  array  of  equally  spaced  elements. 

The  characteristics  of  the  various  tracking  techniques  were  summarized 
after  a  study  of  the  literature.*  The  following  section  presents  comparative  dis¬ 
cussions  of  the  techniques  studied,  and  summaries  of  their  major  features. 

The  main  basis  for  comparison  in  the  study  was  an  additive  system  which 
forms  a  narrow  beam  by  the  summation  in  phase  of  all  the  outputs  of  the  array  ele¬ 
ments.  Of  particular  interest  are  sum  and  difference  systems  such  as  those  used 
in  monopulse  tracking  and  the  "split-beam"  multiplicative  arrays  often  used  in 
passive  systems. 

Of  special  concern  in  the  study  were  the  directional  pattern,  noise  factor, 
and  tracking  sensitivity.  The  directional  pattern  defines  the  one-way  array  gain 
when  evaluated  along  the  bores ight.  and  the  main  lobe  size  as  measured  between 
the  first  set  of  nulls.  In  general,  each  system  will  provide  an  even-order  pattern 
with  a  maximum  at  bores  ight  lor  detection  and  an  odd-order  pattern  with  a  null  at 
bores  ight  for  use  in  tracking  the  detected  target. 

One  other  form  of  the  phase  for  the  ft11'  element  may  l>e  useful.  If  the 
difference  of  direction  cosines  are  replaced  by  the  rr.  r  ,  and  rz  terms  and  the 
location  of  the  /,,th  elements  is  left  in  rectangular  coordinates,  then  6/{  is  given  by 


k 


*  Specific  references  are  riled  in  the  discussion;  other  pertinent  literature  is  ,ils.>  included 
III  I  he  References  section,  pu*»e  til. 


Additive  Linear  Array 


The  general  configuration  of  an  additive  line  array,  which  will  lx?  used  as 
the  basis  of  comparison  between  systems,  is  shown  in  figure  14.  The  overall 
length  of  the  array  is  D  with  n  elements  spaced  d  apart.  With  an  input  signal  of 
.\  cos  the  summing  of  the  n  equally  weighted  transducers  outputs  gives  the 
following" 


COS  ajot 


The  symbol  shall  lx?  used  for  the  output  of  each  beam-forming  channel  and  A  for 
the  output  of  the  nulling  channels  even  though  they  are  most  often  associated  with 
the  sum  and  difference  channels  of  monopulse  systems.  The  context  will  make 
clear  the  system  being  considered. 


Figure  14.  Block  diagram  of  additive  line  array. 


I’sintf  th(*  variable  p  "  sin  0  and  noting  that  for  closely  spac'd  iMrimTii 
tht‘  situ*  function  in  the  denominator  may  \n:  replaced  by  its  argument  tfiws 


V  An  (^LIU>)  COS  t!,J 

\  up  ! 


Tilt'  introduction  of  time  delays  between  the  elements  and  summing  network  pro¬ 
vides  tor  steering  the  main  lx*am.  Since  the  only  change  in  the  expression  given 
is  that  p  is  replaced  by  pl  where  />'  rrti  ,v(sin  n  sin  and  I)  is  the  steering 
angle  introduced,  the  effects  of  beam  steering  will  not  be  included. 

When  a  lxuim-forming  system  is  used  to  estimate  the  angular  loeatiou  of 
the  target  as  well  as  to  detect  the  target,  the  best  estimate  of  target  liearing  is  the 
boresight  of  the  lx.’ am  with  a  standard  deviation  of  approximately  the  half-power 
lieam  width  (half  angle).  Thus  very  little  can  lie  said  about  the  liearing  except  that 
the  target  is  contained  within  the  lieant  (assuming  that  the  side  lobes  have  Ix'cn 
redact'd  to  a  negligible  level). 

In  radar  and  some  types  of  sonar  systems.  ;>  moving  antenna  lolie  is  used 
to  improve  the  angular  accuracy.  Many  types  of  lolie-switehing  and  'ocam-scanning 
tracking  systems  have  lieen  develop'd.  ail  of  which  have  the  common  disadvantage 
of  requiring  several  samples  of  target  data  to  give  an  estimate  of  the'  angle  of 
arrival.  This  is  undesirable  in  systems  where  the  target  returns  undergo  fluctu¬ 
ations  Ix'tween  sample  points.  In  the  case  of  sonar  these  fluctuations  art'  due  10 
target  motion  Ix'tween  samples,  changes  in  the  propagation  path,  platform  motion. 


Systems  which  obtain  the  information  necessary  to  make  an  estimate  of 
the  angle  of  ai rival  on  the  basis  of  a  single  sample  or  pulse  are  classed  as  simul¬ 
taneous-lolling,  or  monopulse,  systems. 

Figure  lr>  shows  one'  way  of  processing  the  outputs  of  the  linear  array  to 
obtain  a  better  estimate  of  the  target  bearing.  The  array  is  divided  into  two  equal 
sulsarrays  where  outputs  arc  summed  and  differenced.  The  sum  channel  will  pro¬ 
vide  a  Ix'am  at  boresight  for  use  in  detection  and  the  difference  channel  will  give 
a  null  at  boresight  to  Ik ■  used  in  tracking. 

Tlx*  outputs  of  each  ef  the  still-arrays  arc  given  by 


where 


<£  =  sm  &  =  nP 


The  sum  of  the  two  sub-arrays  is 


1  =  An 


sin  np 
np 


COS  (Oof 


which,  as  would  be  expected,  is  the  response  previously  stated  for  the  summed 
output  for  the  full  array.  The  difference  channel  output  is  given  by 


The  difference  function,  A,  is  an  odd  function  with  a  null  at  p  =  0  which  makes  it 
a  useful  function  for  tracking  a  target  in  angle.  Since  the  accuracy  with  which  the 
null  point  can  be  determined  is  limited  by  the  noise  present  in  the  difference  chan¬ 
nel  and  is  proportional  to  the  slope  of  the  function  at  the  null,  a  measure  of  the 
system  sensitivity  is  the  ratio  of  the  slope  of  the  direction  pattern  at  p  =  0  to  the 
rms  noise  level  in  the  difference  channel  for  some  normalized  input  signal.  The 
slope  of  the  split-beam  additive  system  difference  channel  envelope  is 


limit  ef  A  _  lim  .  .  P"(sinf  )M)~hf )' 

p  -0  dp  p-»0  p2 


A'  = 


n  2  A 
2 


Thus  the  sum  channel  may  be  used  to  detect  the  signal  and  the  difference  channel 
used  to  track  the  target  within  the  beam  width  of  the  sum  channel.  If  it  is  desired 
to  estimate  the  position  of  the  target  relative  to  the  boresight  rather  than  maintain 
the  target  at  the  null,  it  is  necessary  to  eliminate  the  dependence  on  the  signal 
amplitude  in  the  difference  channel.  It  was  shown  that  for  the  additive  linear  array 
the  difference  signal  and  the  slope  of  the  difference  signal  with  respect  to  angle 
of  arrival  are  both  functions  of  the  input  signal  amplitude.  To  eliminate  this  de¬ 
pendence  on  signal  amplitude  the  sum-channel  signal  may  be  used  to  control  an 
AGC  loop  in  the  difference  channel.  A  system  of  this  type  is  known  as  a  phase¬ 
sensing  monopulse  system. 


Monopulse  Tracking  Systems 

The  two  basic  types  of  monopulse  systems  are  classed  as  amplitude  • 
sensing  and  phase  sensing.  They  are  distinguished  by  the  type  of  antenna  sys¬ 
tem  used  to  produce  the  sum  and  difference  signals.  After  the  sum  and  difference 
signals  are  produced,  their  processing  may  be  the  same  in  either  system.  A  sum¬ 
mary  of  each  type  will  be  given  here;  a  detailed  development  and  analysis  of  each 
may  be  found  in  references  9,  10,  and  11. 

Two  of  the  requirements  of  a  true  monopulse  system  are  that  an  odd  func¬ 
tion  ot  the  angle  of  arrival  be  formed  and  that  this  function  be  normalized  so  that 
it  is  not  a  function  of  the  absolute  value  of  the  received  signal.  By  using  only 
the  ratio  of  the  signal  amplitudes  the  system  is  made  insensitive  to  target-strength 
fluctuations  due  to  variations  in  the  propagation  media  or  target  scintillation. 

In  the  amplitude-sensing  system  (fig.  IK)  two  beams  are  formed  with  the 
axis  symmetrically  displaced  from  the  desired  boresight  by  a  squint  angle  <)s.  The 
phase  centers  of  each  lobe  are  made  to  coincide  and  the  amplitude  patterns  ;iro 
mirror  images  of  one  another  about  the  boresight.  When  using  an  array  of  trans¬ 
ducers  the  pair  of  overlapping  lobes  may  be  produced  in  the  beam  former  with  a 
power  splitter  and  two  fixed  delays  for  each  element  prior  to  summing,  as  shown 
in  figure  16.  Thus,  the  beam-steering  delays  are  not  affected  by  the  use  of 
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ONE  INPUT  TO 
BEAM-FORNER  A 


ONE  INPUT  TO 
BEAM-FORNER  B 


OUTPUT  FROM 
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PHASE 

REVERSAL 
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Figure  16.  Amplitude-sensing  monopulse 


monopulse  tracking.  For  the  case  of  uniform  weighting  of  the  elements,  the  outputs 
of  the  adjacent  beams  are  given  by 

E,  =  D  (p,)  .-A  co.Sf „0t 

Ej  -  D  (;>i )  A  cos  w0 / 

where 

sinnp,-  sinnp, 

— — 1  /=  1.2 

;;sinp(  np, 

Pi  =  ~  (sint?- sin  fls.) 

A 

P2  =  (sin  t?-  sin  ds.) 

A 

The  sum  and  difference  signals  are 

i  =  E,  ■  E,  =  .A  [Dtp,)  •  Dtp,)]  cos w0/ 

A  ■-  E,  E,  r =  .A [Dtp,)  Dtp,)] cos w0/ 

The  normalized  difference  signal  is 

A  _  [Dtp,)  -  D(Pi>] 

A  [Dtp,)  -  Dtp,)] 

f'igure  17  is  a  plot  of  the  sum  and  difference  patterns  vs  p  for  a  cosine 
illuminated  array.  It  is  noted  that  one  important  parameter  of  the  amplitude  mono¬ 
pulse  system  is  the  squint  angle.  0s.  The  squint  angle  is  usually  chosen  to  give 
maximum  linearity  of  the  angle  output,  or  for  maximum  sensitivity  on  the  boresight. 
for  a  given  illumination  function.  The  Kerr-Murdock  condition  11  shows  that  if  the 
illumination  function  is  a  positive  power  of  a  half-cycle  cosine,  a  squint  angle  of 
np  =  77/2  gives  an  output  which  is  linear  with  angle  off  boresight.  For  the  case 
where  it  is  desired  to  maximize  the  boresight  sensitivity  the  choice  is  not  as 
simple.  The  sensitivity  on  boresight  of  the  normalized  difference  signal  is  given 
by  A 10)/ -(0).  Figures  18  and  10  show  a  plot  of  this  function  vs  squint  angle  for 
uniformly  illuminated  and  half-cycle  cosine  illuminated  apertures.  In  both  cases 
the  sensitivity  is  a  monotonically  increasing  function  of  squint  angle  and  thus  has 
no  maximum  within  the  sum  beam  width.  In  radar  and  sonar  applications  where  the 
target  is  illuminated  by  the  sum  beam,  there  is  logic  in  the  selection  of  a  squint 
angle  which  gives  a  maximum  sum  slope  product.  As  is  seen  in  the  figures  these 
maxima  do  exist  and  are  found  to  be  very  close  to  the  3-dB  power  point  independent 
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of  the  type  of  illumination  function  chosen.  For  uniform  illumination  the  maximum 
is  at  np  =  1.30  which  is  2.6  dB  down  on  the  sum  pattern.  The  3-dB  point  is  np  = 
1.39  and  the  6-dB  point  is  at  1.894.  It  was  found  that  the  sum  slope  product  at  the 
6-dB  point  is  about  0.8  of  the  sum  slope  product  at  the  3-dB  point,  or  a  loss  of  less 
than  1  dB. 


-An  -2n  -n  0  n  2n  4?r 


np 

Figure  17.  Sum  and  difference  functions  for  a 
cosine-weighted  additive  array  with  0S  = 


Figure  18.  Functions  of  squint  angle  for  uniform  illumination. 


Figure  1H.  Functions  of  squint  angle  for  cosine  illumination. 


Figure  20  is  a  block  diagram  of  a  processor  for  the  sum  and  difference 
channels.  Both  channels  are  heterodyned  to  an  i-f  where  they  are  amplified.  The 
gain  of  the  i-f  amplifiers  is  controlled  by  the  signal  in  the  sum  channel  through  the 
instantaneous  automatic  gain  control  (1AGC).  Both  channels  are  then  applied  to 
a  phase  comparator  or  phase-sensitivity  amplitude  detector  which  provides  an  out¬ 
put  proportional  to  the  difference  over  the  sum  and  with  the  sign  indicating  sense. 


A  *  p  > 
<p] 


Figure  20.  Example  of  monopulse  processor  for  sum  and  difference  channels. 
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The  rms  error  in  the  estimate  of  angle  of  arrival  using  the  monopulse  will  be  many 
times  less  than  the  beam  width  of  the  sum  channel  for  large  signal-to-noise  ratios. 
An  estimate  of  that  error  is  given  by 


<70  = 


where 


Oo  is  the  half-power  beam  width 
c 

—  is  the  i-f'  signal-to-noise  power  ratio 
A 

K  is  the  slope  of  the  normalized  difference  signal  in  rms  volts  per 
beam  width. 

For  the  amplitude  monopulse  system  with  the  squint  angle  equal  to  the 
3-dB  beam  width,  the  slope  is  approximately  1.5  with  uniform  weighting.  Thus, 
the  normalized  rms  angular  variance  is 


N 

2K3S 


°-222(f) 


The  phase-sensing  monopulse  system  also  forms  sum  and  difference  sig¬ 
nals  which  may  be  processed  in  the  same  manner  as  for  the  amplitude-sensing 
system  previously  described.  The  difference  is  in  the  way  the  sum  and  difference 
signals  are  obtained.  Figure  15  shows  an  additive  array  which  forms  a  sum  and 
difference  signal  by  splitting  the  array  into  two  sub-arrays.  Since  the  arrays  are 
the  same  size  and  both  have  parallel  boresights,  the  far-field  amplitude  patterns 
will  differ  very  little  and  only  the  phase  of  the  sub-array  outputs  will  differ,  with 
the  phase  being  a  function  of  the  angle  of  arrival.  The  amplitude  of  the  i-f  signal 
in  the  sum  and  difference  channels  was  shown  to  be 


^sirwrp 

np 


For  a  fixed-length  array,  greater  displacement  of  phase  centers  can  be  achieved 
by  using  only  part  of  the  array  located  near  the  extremes  with  an  associated 


increase  in  sensitivity  to  angle  of  arrival.  However  the  maximum  amplitude  of  the 
sum  channel  will  be  reduced  in  proportion  to  the  percent  of  the  array  not  used,  and 
Rhodes  11  has  shown  that  the  maximum  sum  slope  product  is  achieved  when  the  full 
array  is  used  with  a  phase  center  separation  equal  to  half  the  array  length.  The 
normalized  slope,  h.  used  in  estimating  the  rms  angular  error  for  tho  phase-sensing 
monopulse  is  approximately  1.4  for  uniform  weighting  of  the  elements.  Thus  the 
normalized  variance  is 


Some  comments  on  the  effects  of  phase  and  amplitude  errors  in  the  two  types  of 
systems  may  be  helpful.  In  the  amplitude-sensing  system,  any  voltage  unbalance 
prior  to  forming  the  sum  and  difference  signals  will  cause  a  shift  in  the  boresight: 
in  the  phase  system  any  phase  shift  prior  to  forming  sum  and  difference  will  cause 
a  boresight  shift.  A  shift  in  phase  with  the  amplitude-sensing  system  will  reduce 
the  null  depth,  w’hile  with  the  phase-sensing  system  any  voltage  unbalance  will 
cause  a  reduction  in  null  depth.  Phase  shifts  and  voltage  unbalances  after  the 
formation  of  the  sum  and  difference  signals  will  not  affect  the  position  or  null 
depth,  but  will  influence  the  sensitivity  of  the  detector.’ 


Multiplicative  Arrays 

Multiplicative  or  cross-correlation  receivers  have  been  used  to  great  ad¬ 
vantage  in  many  Helds  as  a  means  of  reducing  the  effective  beam  width  of  a  system 
and  thereby  improving  the  angular  resolution  (ability  to  distinguish  two  adjacent 
targets).  Most  of  these  applications  have  involved  (1)  very  long  integration  times, 
where  the  signal-to-noise  ratio  may  be  made  as  large  as  desired,  and  (2)  passive 
detection  techniques.  The  performance  of  the  multiplicative  array  will  be  com¬ 
pared  here  with  that  of  an  additive  array,  for  application  to  an  active  sonar  or 
radar  system. 

In  general,  a  line  array  of  n  elements  may  be  combined  in  many  ways  by 
using  only  multiplication  and  averaging.  Many  interesting  and  unique  directional 
pattern?  may  be  synthesized  by  this  method,  but  all  have  the  major  disadvantage 
of  having  an  output  signal-to-noise  ratio  which  is  generally  poorer  than  that  of  the 
smallest  sub-array  used.  In  addition,  the  response  to  multiple  targets  cannot  in 
general  lx?  related  to  the  position  and  relative  amplitude  of  each  target,  because  of 
the  cross  terms  which  arise  in  the  multiplication.  Reference  12  has  a  discussion 
for  the  response  for  two  targets  when  a  single  multiplication  is  used.  If  only  a 
single  multiplication  is  allowed,  an  array  of  n  elements  may  be  divided  into  two 
sub-arrays  of  n,  and  n2  elements.  Each  sub-array  is  operated  as  an  additive  line 
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The  phase  difference  between  the  outputs  of  each  array  is  a  function  of  the  dis¬ 
tance  between  phase  centers,  s.  and  the  angle  of  arrival  of  the  plane  wave  relative 
to  broadside  (no  steering  delays  applied).  Thus 


Xi  =  - 


h 

9 


2  nd 

A 


sin  0 


-Y>  -  X'j 


nd 

2  \  A 


np 


Tlie  output  after  multiplying  and  averaging  as  a  function  of  p  is  proportional  to 


cos  np 


or  since  nt  ■  n2  =?  n  and  sin  p-p  for  small  element  spacing 


Note  that  the  location  of  the  first  null  is  determined  by  the  cos  np  term  which  is 
not  a  function  of  n,  or  n}.  The  selection  of  n,  and  n2  docs  affect  the  shape  and 
side-lobe  levels  of  the  directional  pattern  as  well  as  the  signal-to-noise  perform¬ 
ance.  (Reference  12  discusses  the  effect  on  side-lobe  level.)  The  directional 
pattern  of  an  additive  array  of  n  elements  was  shown  to  lx?  proportional  to 


Dn(p) 


sin  np 
np 


Note  that  the  first  null  for  the  additive  array  occurs  at  np  -  -  while  the  first  null 
for  the  multiplicative  array  is  located  at  np  -  -  2.  Thus  the  angular  beam  width 
to  the  first  null  of  the  multiplicative  array  is  half  that  of  the  additive  array  of  the 
same  length. 

In  considering  the  effect  of  uncorrelated  background  noise  on  the  reception 
of  a  single  coherent  signal.  Welsbv  and  Tucker  13  have  shown  that  the  rms  signal- 
to-noise  ratio.  R  (defined  as  the  ratio  of  the  de  terms  to  the  ac  terms  out  of  the 
multiplier),  is  given  by 

D  D 

R  = - ^ - - 

[Ml  .R.’-R^IS 
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where  ft,  and  ft,  are  the  mis  signal-to-noise  ratios  at  the  two  inputs  to  the  multi¬ 
plier.  If  ft o  denotes  the  mis  signal-to-noise  ratio  at  each  element  and  all  element: 
have  equal  weighting,  then  R ,  s/iTl  K„  and  R ,  v/.<r;  ft,,  and 


R 


\Zn7n,  ftp2 _ 

[S  ( 1  •  n,  ft02  -  n:  Ro2>llj 


If  ft,  and  R2  are  large  eompared  to  i,  then 


or  since  n,  ■  n,  n 


R  -  Rn)fn  / - '=■ 

~  ' 

Vh 

But\/n  R o  is  iust  the  signal-to-noise  ratio  for  an  additive  array  of  a  elements;  thus 
the  noise  figure  of  the  multiplicative  array  is  given  by 


Figure  22  is  a  plot  of  the  noise  figure  as  a  function  of  the  ratio  of  a,  a,  assuming 
ft,2  *  ft,-  -  l.  The  minimum  noise  figure  is  equal  to  \fl  when  a,  a,  1.  Thus, 
the  multiplicative  array  has  a  signal-to-noise  degradation  of  at  least  :l  dB  com¬ 
pared  to  an  additive  array  of  equal  length.  As  the  ratio  of  a,  a.  Iiecomes  large  the 
noise  figure  approaches  (a  2>V  Thus  when  a,  1  the  noise  figure  Iiecomes  pro- 
giessively  worse  as  the  numher  of  elements  is  increased.  This  is  because  the 
signal-plus-noise  for  the  single  element  is  always  one  of  the  inputs  to  the  multi¬ 
plier  and  output  signal-to-noise  ratio  cannot  he  greater  than  that  of  the  single 
element. 

Setting  a,  a,  a  2  to  maximize  the  signal-to-noise  ratio,  the  output  of 
each  sub-array  due  to  a  signal  \  eos  insistent  at  an  angle  H  (fig.  2'i)  is  given 


Figure  24  is  a  plot  of  the  array  response  vs  the  parameter  p  for  the  split- 
beam  multiplicative  array  with  nl/nl  =  1  and  n,/n3  =  8  compared  with  that  of  a 
uniform  additive  array. 


Figure  24.  1  pattern  for  multiplicative  arrays. 

For  hearing  determination  using  a  split-beam  multiplicative  system,  a 
phase  shift  of  .t  2  radians  is  introduced  in  one  channel  before  multiplication  to 
provide  a  signal  which  is  an  odd  function  of  angle  of  arrival  and  has  a  null  at 
boresight: 


This  function  is  plotted  in  figure  25  with  An  =  1.  The  slope  of  this  function  at 
boresight  (p  =  0)  is  twice  that  of  the  split-beam  additive  array  null  function.  Since 
the  signal-to-noise  voltage  ratio  is  down  by  the  \[2 .  the  multiplicative  split-beam 
system  is  the  \  2  more  sensitive  to  arrival  angle  near  the  boresight  when  using 
normalized  sum-and-difference  systems. 
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Figure  25.  A  pattern  for  split-beam 
multiplicative  array. 


Interferometers 

For  the  purpose  of  completeness,  the  characteristics  of  the  two-element 
additive  and  multiplicative  interferometers  will  be  given.  Figures  26  and  27  show 
a  two-element  interferometer  separated  by  a  distance  D  =  nd.  The  output  of  each 
element  as  a  function  of  the  angle  of  arrival  is  given  by 


Ei  =  A  cos  W  +  y) 
P 
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For  the  additive  interferometer  (fig.  26)  the  sum  channel  output  is  given  by 
I  =  E,  +  E,  =  .4  2  (cos  coot)  (cos  x>] 

1  =  2 A  cosnp^  cos  coDt 

and  the  difference  channel  by 

A  =  2 A  (“~y2(sin<y0f)(sin-x)] 

A  =  -  2 A  ^2£I?  •  sinnpj  sinw0< 

Thus  the  additive  interferometer  has  a  main  beam  width  half  that  of  a  full 
additive  array  whose  length  is  equal  to  the  interferometer  spacing.  This  is  because 
the  phase  centers  are  approximately  twice  as  far  apart.  However,  the  signal-to- 
noise  performance  is  down  by  \[n/ 2  because  of  the  use  of  fewer  elements,  and 
grating  lobes  of  almost  equal  gain  as  the  main  lobe  will  appear  at  np  =  kn,  k  =  1, 

2,  3,  ••  •  giving  angular  ambiguities. 

For  the  multiplicative  interferometer  the  beam-forming  channel  output 
(fig.  27)  will  be 

1  =  EtE}  =  A2(^py(V2  cos 2*) 

A’/sinpV  o 
=  t  (  — - )  cos  2  np 
2  \  p  / 

The  null  channel  will  give 


Thus  the  multiplicative  interferometer  reduces  beam  width  by  a  factor  of  2  and 
gives  twice  as  many  grating  lobes  as  the  additive  interferometer.  The  power  sig- 
nal-to-noise  ratio  for  the  multiplicative  system  is  reduced  by  3  dB  over  the  ad¬ 
ditive  system,  as  expected. 
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Comparison  cf  Systems 


A  summary  of  the  characteristics  of  the  various  systems  is  contained  in 
tables  2  and  3.  Table  2  gives  the  directional  patterns  for  the  Ixuun-forming  sys¬ 
tems  where  the  amplitude  of  the  additive  array  has  been  normalized  to  1  at  the 
j)eak  response  and  all  other  systems  are  referenced  to  the  additive  system.  The 
location  of  the  first  null  is  used  as  an  indication  of  the  main  lolx>  width  and  the 
noise  factor  is  defined  as  the  ratio  of  the  signal-to-noise  for  the  additive  system 
over  that  of  the  system  considered.  Thus  an  increase  in  the  noise  factor  repre¬ 
sents  a  poorer  signal-to-noise  performance  for  a  coherent  signal  and  for  uniformly 
distributed,  uncorrelated  noise. 


TABLE  2.  CHARACTERISTICS  OF  BEAM-FORMING  SYSTEM 


Technique 

Directional  Pattern 

Location  of 
1st  Null 

Noise 

Factor 

(ratio) 

Noise 

Factor 

(dB) 

Uniformly  illuminated 

additive  array 

' 

sin  np 
rip 

TT 

P  =- 
n 

1 

(ref.) 

0 

Cosine  taper  additive 
array 

sin  |/ip  -  ^  sin  (rip 

j 

D  =  3»r 

1.11 

o  kn 

2  cos  — 

n 

h,  0,  ±  i,  ±2,  •  •  ±n  2 

j 

[np 

2/  1 

1  77  \ 

np  +  2j 

'In 

Spl  it-beam 
multiplicative  array 

1  / 

fs\r,  as\ 

2 

1 

..  TT 

| 

v  2 

S.\  •  i 

— 

— 

|  cos  np 

P  -  7T 

3.0 

«  \ 

V  t/ 

In 

motiopulse  system.  The  S  .V  ratio  is  important  in  locating  the  null  since  the  ii-.msc 
in  the  null  channel  is  the  factor  which  produces  an  error  in  estimating'  the  angle. 

A  comparative  ratio  similar  to  that  defining'  the  noise  figure  is  given  by 


r 


where  \  is  the  null  channel  directional  pattern  and  A  is  the  beam-forming  channel 
directional  pattern.  The  numerator  is  related  to  the  additive  array  which  is  used 


TABLE  T  CHARACTERISTICS  OF  NULL  SYSTEM 


System 

ji  (rate) 

Utratio) 

U(dB) 

Uniformly  illuminated 
difference  system 

mj 

1 
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0 

np 
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k  0,  ±1,  ±2.  •  •  ■  tn  2 
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interferometer 
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SUMMARY  AND  CONCLUSIONS 


A  comparison  of’ the  techniques  for  processing  the  information  from  a  line 
array  of  elements  has  shown  that  it  is  possible  to  obtain  improvement  in  one  sys¬ 
tem  characteristic  (i.e.,  angular  resolution)  at  the  expense  of  another  <i.e..  signal- 
to-noise).  The  criterion  for  selecting  any  particular  technique  is  the  importance 
of  each  system  parameter  to  overall  system  performance.  If.  as  is  often  the  case, 
the  maximization  of  signal-to-noise  ratio  is  most  important,  it  was  shown  that  a 
uniformly  illuminated  additive  array  provides  the  best  signal-to-noise  ratio  for  a 
given  array  size  and  for  uncorrelated  uniformly  distributed  noise  background.  Any 
reduction  in  signal-to-noise  ratio  will  result  in  a  reduced  probability  of  detection. 
Thus,  for  detection  in  a  noise  background,  the  additive  array  is  lx>st.  If,  on  the 
other  hand,  some  sacrifice  in  detection  probability  can  lie  tolerated  in  exchange 
for  improved  angular  resolution,  a  split-beam  multiplicative  system  provides  a 
main-lobe  beam  width  half  as  wide  as  that  of  an  additive  array  of  the  same  length 
with  a  ioss  of  3  dC  in  signal-to-noise  powei  ratio. 

The  use  of  only  the  two  e.id  elements  of  tiie  array  in  an  additive  inter¬ 
ferometer  provides  the  same  angular  resolution  that  is  possible  with  the  split-beam 
multiplicative  array,  but  at  a  great  ioss  in  signal-to-noise  and  with  the  addition  of 
ambiguities  in  angle  due  to  the  grating  lobes.  Tne  use  of  the  same  two  end  ele¬ 
ments  in  a  multiplicative  interferometer  provides  improvement  in  angular  resolution 
Coy  a  factor  of  2)  over  the  addit  ive  interferometer,  in  exchange  for  a  3-d B  loss  in 
signal-to-noise  ratio  and  a  doubling  of  the  number  of  grating  lobes. 

The  multiplicative  system  produces  cross  terms  which,  in  the  ease  of 
multiple  targets,  may  result  in  outputs  tuat  have  no  correlation  with  the  target 
positions.  To  prevent  these  undesired  outputs,  range  gating  within  each  beam  will 
be  necessary  to  provide  resolution  in  range. 

To  improve  the  accuracy  of  the  estimate  of  bearing,  a  “nulling”  system 
may  lx?  used  with  any  of  the  beam  forming  systems  discussed.  In  general  the 
nulling  systems  are  the  same  as  the  beam-forming  systems,  except  that  in  the  ad¬ 
ditive  systems  a  phase  shift,  of  n  radians  is  introduced  in  one  arm  prior  to  summing 
and  in  the  multiplicative  systems  a  shift  of  t ' 2  radians  is  made  prior  to  multi¬ 
plying.  In  comparing  the  nulling  systems,  two  parameters  are  important.  The 
slope  of  the  function  near  the  null  is  related  both  to  the  sensitivity  to  displace¬ 
ments  from  ;he  null  and  to  the  rms  value  of  the  noise  present  in  the  null  channel 
tor  some  standard  input,  since  at  the  null  the  noise  level  produces  the  error  in  the 
angular  estimate.  In  table  3  the  directional  patterns  of  the  nulling  systems  as¬ 
sociated  with  the  Ixiam- forming  systems  of  table  2  are  given.  The  ratio  ft  is  the 
slope  of  the  directional  function  at  the  null  normalized  by  the  peak  signal  from  the 
Ixunn- forming  function.  These  have  Iteen  compared  with  the  uniformly  illuminated 
additive  array  again  and  the  smaller  the  ratio  the  greater  the  slope.  The  furictio" 

1  is  equal  to  ft  times  the  ratio  of  the  signal-to-uoi.su  ratios  (the  noise  factor 


given  in  table  '*).  Table  II  shows  that  the  sine  taper  provides  slightly  (-1.2  (IB) 
better  angular  accuracy  despite  the  slightly  (0.92  (IB)  poor  signal-tonoise  per¬ 
formance .  The  split-beam  multiplicative  array  provides  a  factor  of\  2  improve¬ 
ment  in  angular  accuracy  and  a  factor  of  2  improvement  in  angular  resolution,  in 
exchange  for  the  3-dB  loss  in  signai-to-noise  ratio.  All  the  nulling  systems  as¬ 
sume  that  only  a  single  target  is  present  within  a  resolution  cell. 

It  appears  tiiat  the  small  improvement  in  angular  accuracy  and  resolution 
available  with  the  split-beam  multiplicative  systems  are  not  worth  the  price  of  a 
•  5-dB  loss  in  signal-to-noise  in  an  ASVV  situation  where  detection  is  of  prime  im¬ 
portance. 
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